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Abstract 

In this work the qualitative analysis of statics and dynamics of de- 
fects and textures in liquid crystals is performed with help of meanders 
and train tracks. It is argued that similar analysis can be applied to 
2+1 gravity.More rigorous mathematical justifications are presented in 
the companion paper (Part II) on quadratic differentials and measured 
foliations. Meanders were recently introduced by V.Arnold (Siberian 
J. of Mathematics 29,36 (1988)) and are used originally in the combi- 
natorial problem of finding the number of distinct ways given curve can 
intersect another curve in prescribed number of points fixed along this 
auxiliary curve. Train tracks were introduced by W.Thurston (Geome- 
try and Topology of 3-Manifolds, Princeton U. Lecture Notes, 1979) in 
connection with description of homeomorphisms of two dimensional 
surf aces. Train tracks alone are sufficient for the description of stat- 
ics and dynamics of liquid crystals and gravity .Using train tracks 
the master equation is obtained which could be used alternatively to 
the Wheeler-DeWitt equation for 2+1 gravity.Since solution of this 
equation is possible but requires large scale numerical work, in this 
paper we resort to the approximation of train tracks by the meanditic 
labyrinths. This then allows us to analyse possible phases (and phase 
transitions) of gravity and liquid crystals using Peierls-like arguments 
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1 Introduction 



1.1 Motivations and background 

The role of topology in solving physical problems is steadily increasing [1,2 
].In some instances it has become difficult to decide wether physical argu- 
ments are helping to solve topological problems or topological arguments 
are helping to solve physical problems [2]. Here we would like to provide yet 
another example of this sort. 

In 3+1 dimensions the laws of Newton's gravity and Coulombic electro- 
static look suspiciously similar:both are being described by the Poisson-type 
equation 

W = --p, (i.i) 

e 

where e is the dielectric constant in the case of electrostatics while e~ 1 =G 
is the gravitational constant in the case of gravity. The density of charges 
p(r)can be both positive and negative for electrostatics while only positive 
in the case of gravity. The solution of Eq.(l.l) for the potential ip(r) is given 
by 

<p(r) = - [ d\- - P(r ° = (1.2) 

where r={x,y,z}and V is the volume which encloses the charges.lt is ex- 
pected, that if the charges are confined within the domain V,then <p(r) ^0 
when |r| — > oo. Indeed, let p(r) = g5(r),then we obtain trivialy 

<p(r) = l± (1.3) 

where q is the magnitude of charge placed at the orign. Solution given by 
Eq.(1.2) is valid only in 3 dimensions,however.In two dimensions it should 
be replaced by 

VOO = \ I d 2 r'p (r') ln( 1 ==) (1.4) 

e J y/( x - x' f + (y- y') 2 
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(with volume V being replaced by the area A) so that, instead of (1.3) ,we 
obtain now 

ip(r) = -^lnlrl . (1.5) 

£ 

This time, for |r| — > oo we obtain </?(r) — > ±00 (depending upon the sign 
of q).The result (1.5) can be obtained directly from Eq. (1.2) for a special 
case of charges uniformly distributed along the infinitely long and infinitely 
thin rod placed perpendicular to the x-y plane, i.e. 

-- \n(J(x - x') 2 + (y- y') 2 ) = lim (- f dz\ - - - — - --m2L) 

2 vvv J vy yjJ l^oo K 2J V (x-x') 2 + (y-y') 2 + (z - z') 2 ' 

-L V 

(1.6) 

Instead of a rod we can imagine a particle "living" in 2+1 dimensions(i.e. particle 
which evolves in time while moving in x-y plane). 

The difference of behaviors of potentials at infinity in two and three di- 
mensions leads to the profound difference in the underlying physics. More- 
over, this difference forbids smooth dimensional continuation of the expression 
for the potential, given in d-dimensions by 

Mr) = I -U (1.7) 
l r l 

x a — 1 

based on the well known limiting procedure ,e.g. lnx= lim . Let 

a->0 a 

us explain why this is so. Since the potential (p is an auxiliary quantity in 
both gravity and electrostatics,we shall use instead the forces F obtained in 
a usual way through 

F = -Vv? • (1.8) 

The forces create the vector fields and, in the case of 2 dimensions, we 
have to consider the vector fields on surfaces. In particular, let us consider 
the plane R 2 ,the sphere S 2 = R 2 U {00} and the disc D 2 of some radius 
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R.Evidently,physically all three possibilities are almost indistinguishable:the 
sphere S 2 can be obtained from the plane by one point compactification, the 
plane can be obtained from the disc D 2 by taking the limit oo .Topolog- 
ically,however,the above tree surfaces are completely different: the sphere S 2 
has the Euler characteristic x(S 2 ) = 2, the disc x(D 2 ) = 1 and the plane R 2 
cannot be given any value of \ [3]- For the time being we shall work with S 2 
and D 2 .The choice between these two surfaces is determined by the boundary 
conditions. 

Vector fields on the surface(manifold M) obey the Poincare-Hopf (P-H) 
index theorem [4]: 



X (M)=^J(x 4 ), (1.9) 

i 

where the index I(xj) of isolated singularity is defined as follows. Consider the 
vector field v(x,y)={va;(a;, y),v y (x, y)} around just one of the singularities and 
let U G M be some (say , circular) domain around this singularity. In this 
domain we can construct a unit vector n(x,y) via 



n(x,y) = ^\- . (1.10) 

|v(x,y)| 

This vector provides the Gauss map from M to S 1 .The degree of this mapping 
is I(xj).More explicitly [5,6], 



Kl> 2vr / v 2 

Ci 

I cases of both Newtonian gravity and Coulombic electrostatics the index of 
an isolated singularity is equal to one as can be easily seen by introducing 
the complex variable z = x +y =re ltp and by combining Eqs (1.5), (1.8) and 
(1.11). This means that the minimal number of charges which one can put on 
the sphere is 2 while that on the disc is l.The P-H theorem,Eq.(1.9), does 
not automatically leads to the requirement of electroneutrality: one can 
place two charges of the same sign on S 2 without violating this theorem. The 
complications arise, however, if one would like to place more than 2 charges 
on the sphere. In this case if initially we would have,say,2 pluses,we would be 



(1.11) 
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unable to place the additional charges while if we would have initially "+" 
and "-" charges, then such placement becomes possible if we add charges 
pairvise while keeping the electroneutrality.This becomes possible only be- 
cause the presence of two identical charges creates a saddle (e.g. see Fig. la 
) which has an index -l,Ref.[4]. Hence,in the case of electrostatics the only 
possible singularities are the sources/sinks (Fig. Id )with index 1 and the in- 
duced saddles (Fig.la)with index -1. The existence of multiple charges 
in electrostatics and the electroneutrality are directly related to 
the topology of the underlying manifold and to the emergence of 
the induced saddles. 

Consider now the case of Newton's gravity. Let us place the same 2 
charges (masses) on the sphere and let ,say,one of the charges be " +" so that 
another is "-" and, since this is gravity ,the masses are naturally attracting 
each other. Let us try now to add one additional mass (charge). We imme- 
diately run into problem:since all masses attract each other,there cannot be 
saddles. If this is the case,we are unable to place an additional mass on S 2 
.Hence, in such spherical universe we would have just 2 masses ! The situation 
becomes even more dramatic if, instead of S 2 we would consider D 2 .In this 
case we would not be able to put more than one mass while for the case of an 
annulus(i.e.the disc with a hole)we would not be able to put on it even a single 
mass! since the Euler characteristic of an annulus is 0, Ref.[3]. At the same 
time, we can easily put e.g. 3 charges of the same sign on the annulus in the 
case of electrostatics. These very simple arguments lead us to the conclusion 
that there are some profound differences between the Coulombic electrostat- 
ics and the Newtonian gravity and that the above differences alone are 
sufficient in order to arrive at the correct Einstein formulation of 
gravity.This is going to be demonstrated mainly in the accompanying paper 
[7] while here we provide only the qualitative arguments .For this purpose 
we also need to discuss systems other than electrostatic and gravitational 
(in Newton's sense). These are naturally occurring as defects and textures in 
liquid crystals. 

To study these defects it is helpful to recall some basic facts from the 
qualitative theory of dynamical systems on surfaces [5, 6]. These are also de- 
scribed in terms of flows of the vector fields. In two dimensions we are usually 
dealing with the system of two equations 
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^ = v x {x,y), (1.12) 



dy ( v 
^ = vv(*,y). 



Irrespective to the explicit form of the r.h.s. of Eq.(1.12) ,it is known [5,6] 
that the singularities of the flow could only be of the type depicted on 
Fig. lQ*. These can be easily obtained by linearization of the Eq.(1.12) around Fig.l 
the singularities of the vector field. The whole phase portrait can then be 
built out of these local pictures by gluing together the local pieces in some 
consistent way. Moreover, the topological considerations allow us to restore 
the phase portrait with help of only partial knowlege of the existing sin- 
gularities. This principle is very helpful for our case too but the situation is 
complicated by the fact that, in the case of liquid crystals as well as in 
2+1 gravity , instead of vector fields on surfaces one has to deal with line 
fields. This was demonstrated for liquid crystals in Refs.[ 10,11] and will be 
proven for the case of gravity in Ref.[7] .The case of line fields has been 
also studied in the theory of differential equations on surfaces[5, 6] .In this 
case, instead of the system of Eqs (1.12), we have to consider 

dX dy (1.13) 



Typical singularities of the line field are depicted in Fig. 2 with the values 
of the corresponding indices. For example, in the case of the line fields with Fig. 2 
indices ±~,the differential equation which describes these fields is given (in 
polar r,if coordinates)by[5] 

dr , Lo , . 

— = rtan(n-) . (1.14) 
dip 2 

For n=l we obtain the line field with index -| which is described analyti- 
cally by 

a 



(COS*!)* 



Clearly,the flow directions on some of the figures could be reversed 
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with a being an arbitrary positive constant and ip lying in the following 

sectors 

a)-|<( / 9<|,b)|<v?<7r,c)7r<( / 9<-| !: . 

Similar analysis can be performed for n=-l which produces the index 
+ | as depicted in Fig. 2. Although for the line fields it is impossible 
to introduce consistently the global orientation[4 — 6], it is possible 
, nevertheless, to introduce the indices of these fields [4] by analogy with 
the vector fields. The absence of orientation for the line fields leads to the 
absence of "forces" between the singularities (defects). This property is not 
unusual for the theories of 2+1 gravity where it is impossible to introduce 
the interaction forces between the massive particles[12]. In sections 4,5 of 
Part II (Ref.7) we shall demonstrate that both 2+1 gravity and the textures 
in liquid crystals are described by the line fields. In this part (which we call 
Part I )we shall rely mostly on the intuitive arguments. 

Although it is impossible to introduce the interaction forces between de- 
fects (or masses) it is possible, nevertheless,to talk about the total energy in 
both 2+1 gravity [13] and in the theory of liquid crystals (e.g. see Part II and 
section 5 below ). Moreover, because of this fact, it makes also sense to talk 
about the different phases (orders) at least for the case of defects in liquid 
crystals [14] .Although the connections between the defects and textures in 
liquid crystals and 2+1 gravity were discussed before [15] , here we provide 
completely different treatment of these connections. 

1.2 Organization of the rest of this paper 

In section 2 we consider the role of topology in phase transitions in 2 di- 
mensions. We argue, that in cases of the vector (Coulomb-like) and line fields 
on surfaces topological considerations alone are sufficient for predicting the 
Kosterlitz-Thouless type of phase transition [14] from gas to dipole phase 
while in case of the line fields in liquid crystals ,the existence of hexatic 
phase [16] can be easily established .Clearly , topological considerations alone 
provide only sufficient conditions for existence of ordered phases. The neces- 
sary conditions require us to study the nature of the disordered phase in some 
detail. This is accomplished in sections 3 and 4. In section 3 we provide an 
explicit construction of the disordered phase using some ideas from the qual- 
itative theory of ordinary differential equations. Superposition of these ideas 
with topological arguments produces notions of "labyrinths" and "meanders" 
.Since neither labyrinths nor meanders are widely known in mathematical 
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physics literature ,to avoid repetitions ,we refer the reader directly to section 
3 for precise definitions. Here ,we only would like to mention that the de- 
scription of labyrinths is closely related to the description of the maps of the 
circle [17, 18] which should be more familiar to physically trained reader. The 
notion of a meander has its origin in the work by Arnold [19] who actually 
invented this word. In the simplest case one can think about the meander 
as a representative of the class of curves which intersect another curve in 
the prescribed number of points. Although from the above vague definition 
it may not be clear that the meanders and the labyrinths may have many 
things in common,nevertheless,this is the case as we demonstrate in section 
3. 

In section 4 we discuss yet another way of looking at line fields through 
the notion of "train tracks" . This way of looking at line fields was proposed 
originally by Thurston [20] in connection with his studies of 3-manifolds. We 
argue in this section that the train tracks can be used for description of 
dynamics of textures in liquid crystals and 2+1 gravity.The master equa- 
tion which describes the evolution of the train track could be used as an 
alternative to more familiar Wheeler-DeWitt equation (in the case of grav- 
ity). Since the actual calculations which involve this equation resemble that of 
the Heisenberg-type quantum mechanics, they require a large scale numerical 
work which we try to avoid (at this stage of research) by restricting ourself to 
the approximation of train tracks by meandritic labyrinths. The legitimacy 
of this approximation is discussed in the Appendix. 

Use of the meandritic labyrinth approximation allows us to obtain addi- 
tional quantitative information about the stability of the disordered phase 
and the parameters of the order-disorder phase transition in the system of 
meandritic labyrinths. This is accomplished in section 5. 

All results of Part I are provided without serious proofs (with few ex- 
ceptions). Part II (Ref. [7]) serves to some extent to correct this deficiency by 
providing the mathematical justifications to the emerging picture .Sufficient 
details (and references) are provided in both parts of this work to make this 
presentation self-contained and accessible not only to the experts on gravity 
but to the interested condensed matter researchers as well. 
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2 Role of topology in phase transitions in 2 
dimensions 



We had already demonstrated (in section 1.1.) that the property of elec- 
troneutrality is directly connected with the topology of the underlying man- 
ifold, at least in the case of Coulombic -like systems. Here, we would like 
to demonstrate that ,in addition,the topological considerations alone deter- 
mine the nature of phase transitions in such systems. For the line fields the 
sequence of arguments ,unfortunately,is considerably more complicated and 
will be discussed in the rest of this paper and in Part II. 

In the case of Coulombic-like system, let us assume that under some 
favorable conditions the incoming "+" and "-" have "decided" to stay to- 
gether thus forming a dipole.The vector field in this case is depicted in Fig. 3 
along with a separate source (+) (or sink (-)).The index of the dipole is 2 Fig 3. 
so that P-H theorem is not violated if on the sphere we would have only 2 
charges which have "decided" to form a dipole .Let us place a source/sink in 
the vicinity of such dipole. By direct inspection we can notice that topologi- 
cally such source/sink cannot fit the vector field flow pattern coming from the 
dipole unless it is combined with yet another charge thus forming yet another 
dipole. Formation of the second dipole will cause the formation of 2 additional 
saddles so that P-H theorem is not violated. Evidentlythis procedure can be 
extended to other charges. From here several conclusions could be drawn. 
First, the transition to the dipole phase should be sharp. Second, since in the 
above arguments the density of charges p was not present, it means that 
such transition cannot be considered in traditional thermodynamical sense 
,that is for some critical temperature T c and critical pressure P c there is no 
true critical density p c .These conclusions are in complete accord with the re- 
sults of Hague and Hemmer [21] who had provided more traditional statistical 
mechanics treatment of this type of transition. Obtained results are also in 
complete accord with the results of more sophisticated treatment performed 
by Kosterlitz and Thouless [22]. 

The above simple picture breaks down immediately if we are willing to 
analyze possible transitions in the case of defects which produce the line 
fields, e.g.liquid crystals or 2+1 gravity. Unlike the Coulombic case where the 
vector field picture is rather simple:or we have an independent charges or we 
have dipoles, in the line fields case there are many more possibilities. Since 
the local and global analysis of the line fields had been a subject of intensive 
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research in mathematics for quite some time ,it is impossible to squeeze 
the enormous amount of results accumulated to date in this work. For a 
good summary,please, consult Refs. [23-25 ]. Without even trying to make a 
summary of these results, we shall, nevertheless, select those which ,we feel 
,are of immediate physical relevance. 

Let us begin with the observation that in case of the Coulombic-like fields 
the regrouping of charges had not produced a violation of the P-H theorem. 
Evidentlywe have to require the same in the case of line fields. The moves 
which do not violate P-H theorem are depicted in Fig. 4 and are known in the 
literature as Whitehead moves[23] In addition, it is possible to imagine the Fig.4 
situations depicted in Fig. 5 .Clearlyby looking at Fig.4 we observe the case Fig.5 
of creation of "new " defects with index ±1 from the "old" ones with index 
±|, while looking at Fig. 5 we effectively observe the case of destruction 
of defects. The indices of a) and b) are + \ while the index of c) is zero. 
In principle,other situations are also possible but,for reasons which will be 
explained in section 4 and Part II, we shall not be concerned here with more 
complicated situations. 

Given these local moves, what one can say about the "phase transitions" 
(as compared to the Coulombic -like vector fields) in the system of such 
defects ? For example, can we expect ,by analogy with Coulombic case,that 
the transition(s) is ( are) going to be density-independent? For this purpose 
let us consider the case of a disc x{D 2 ) = l.To be in accord with the P-H 
theorem, Eq.(1.9), we need to have at least four defects with half integer 
indices (since only half integer types of defects are topologically stable [26]) 
:one with index -\ and 3 with the index +| .The resulting "stable phase" 
is depicted in Fig. 6. This does not look like a dipole and,because of the P-H Fig. 6 
theorem,we cannot just form another dipole and so on as in the Coulombic 
case and we cannot use the electroneutrality condition as well. So,what else 
could we imagine? For example, should we have exactly 14 defects ,we could 
have an orderly structure which resembles that depicted in Fig. 7 We use the 
word "periodic" to describe this phase following Thurston's classification of 
surface homeomorphisms [25] as will be explained in more details in section 
4. Fig. 7 

Is it possible to obtain something else ? The answer depends upon how 
many defects are at our disposal. Consider the case of thermodynamic limit 
,i.e.the case when the size R of our disc D 2 is allowed to approach infinity 
.In this situation we can imagine the structure depicted in Fig.8.How can Fig.8 
we check that this structure is actually in agreement with the P-H theorem? 
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This can be accomplished following the original arguments by Hopf 4 (see 
also Ref..27). Let S be closed orientable surface of genus g . Let ao be the 
number of vertices, ai be the number of edges and a 2 be the number of 2-cells 
(e.g.hexagons). Then, according to the Euler theorem,we obtain , 

a - ai + a 2 = 2 - 2g (2.1) 

Looking at Fig. 8 and following the line of arguments made by Hopf,let us 
replace the original hexagonal lattice with the elementary cell depicted in 
Fig.9 Fig.9 

The face of each cell has an index +l,the edge has an index 2(-|) = 
— 1 and each vertex has an index l.Whence,such construction holds for an 
arbitrary g and,therefore,the hexagonal phase is consistent with the Euler 
theorem ,Eq.(1.2) ,and,accordingly,with the P-H theorem,Eq.(1.9).As in the 
case of the Kosterlitz and Thouless transition earlier described ,the existence 
of the above hexatic phase is already known in physics literature,e.g. see 
Ref. [28], where this phase is obtained based on completely different set of 
arguments. In section 4 of Part II we shall provide yet another arguments in 
support of the existence of this phase. 

The questions arise: 

a) is such obtaned hexatic phase unique ? 

b) can such phase occur for an arbitrary concentration of defects (as in 

the 

Coulombic phase case)? 

The answer to the first question is negative as can be easily seen from 
the situation depicted in Fig. 10a) and the corresponding elementary triangle 
given in Fig. 10b). Taking into account that , in view of Fig.5c) ,the effective Fig. 10 
total index of singularities at the edge (aoa2) is zero, we obtain the situation 
which coincides exactly with that described by Hopf [4] .That is we have the 
index of 2-cell the same as a 2 vertex (i.e.+l),the index of the edge the same 
as the index at ai (i.e.-l) and the index of the vertex is the same as the index 
at ao (i.e.+l). Hence, we are again in accord with the Euler theorem, Eq.(2.1). 

Since we have at least two types of hexagonal structures we may think 
about some sort of phase transition between them. In addition, we have not 
actually proved that only two hexagonal structures are possible .A similar 
but much simpler problem of packing of hard hexagons on the triangular 
lattice was considered by Baxter[29]( and, more recently ,by Monasson and 
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Pouliquen) with partial success (since the solution of the hard hexagon model 
is related to the solution of 2d Ising model in the magnetic field which is not 
known in general ). 

Because of some similarities between the hard hexagons and the present 
liquid crystal problem, it is useful to furnish some details since they will 
eventually be used to provide an answer to the second question (about the 
concentration dependence) posed above. To this purpose, in the case of hard 
hexagons ,let us introduce the surface density p via 



(a similar problem which involves the one component plasma was considered 
recently in Refs [30, 31]), where n is the total number of hexagons and N is 
the total number of lattice sites (triangles). More useful quantity ,however,is 
the packing fraction 77 = 3p .To understand the emerging factor of 3 , the 
inspection of Fig.ll is helpful (see also Ref.32 ).By construction, the complete Fig 
packing corresponds to rj = 1 .Looking at Fig .10 ,we easily can conclude 
that, in the case of the packing depicted in fig.lOa),the packing fraction is \ 
.Even in the case of hard hexagons the packing problem is not completely 
under control. That is the description of transition from one mode of partial 
packing to another (e.g. loose r] < 1 versus dense i] = 1) is still lacking [29] 
.At the same time, in the present case we are not even dealing with hard 
hexagons. Our "hexagons" can be easily destroyed as it will be discussed in 
section 6 of Part II. To analyse this more complicated situation, quantities 
other than p may be helpful. Following Ref. [30] ,let us introduce some cut- 
off radius r (the size of the disclination core [14, 26])as well as the surface 
density p via 



P = \ ■ (2.3) 

where A is the surface area and h is the total number of surface defects 
(charges). Next, we can introduce the Wigner -Seitz radius r w _ s via 



7rr„ 



(2.4) 
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Finally,the filling fraction v can now be introduced via 

2 

= ^oP . (2.5) 

By construction, 0< v < 1 so that v can be used (instead of i] )to characterize 
the possible phase transitions. These results will be used in section 5. 

3 Labyrinths and meanders 

The disordered phase of line field defects is no less interesting than the or- 
dered phase. Let us consider the simplest possible case of just four defects in 
D 2 .The "ordered" phase is depicted in Fig. 6 .To get a feeling of the disor- 
dered phase, following Ref.[33], let us consider the "foliation box" JF depicted 
in Fig. 12a). The lines , known in mathematical literature as foliations (mea- Fig. 12 
sured foliations to be exact [10, 23 — 25]) are free of singularities in the box a) 
and contain 2 singularities in the box b)(with account of the results depicted 
in Fig. 5c)). Following Ref. [33] ,we call these singularities "Y" and "thorn" re- 
spectively. Since the total index in the box b) is zero,we can ,evidently,place 
in the line field as many as we wish such Y-thorn doubles without violat- 
ing of the P-H theorem, Eq.(l. 9). This gives us a certain freedom of moving 
such objects around. In particular, if for the moment we would like to forget 
about the "tail " of Y-singularity,then,we can concentrate our attention at 
the foliations in semicircle(s).We shall restore the tail afterwards ( to account 
for the P-H theorem). In the meantime, let us consider collection of 3 foliated 
semicircles as depicted in Fig. 13. Fig. 13 

The foliations for Si are centered at a ,for S2 at b,and for S3 at c respec- 
tively.The ratio of diameters D^/D^ = a is chosen to be some irrational 
number. Let x denote the location of an intersection of a given foliation 
leave with the axis [0,1]. Then, depending upon the direction from which the 
foliation hits the x-axis ,one can distinguish between the following 3 possi- 
bilities: l)if the foliation hits [0,1] by coming from Si(e.g.see point b),then 
let Ri(x) = a — x, < x < a; 2)if the foliation hits [0,1] by coming from 
S 2 ( e.g. see point c ) ,then let R 2 (:r) = 1 — x; 3) if the foliation hits [0,l]by 
coming from S 3 (e.g. see point 2'), then let ~R,3(x) = 1 — (x — a), a < x < 1. 
The meaning of thus introduced functions Ri,R2,and R3 becomes clear if we 
consider their compositions, e.g. [Rio R 2 ](x) and [R 3 o R 2 ](x) . Looking at 
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Fig. 13 ,we can notice ,that in both cases we are dealing with the second 
return map of the closed interval [0.1], e.g. see the paths 1-2-3 (for R^o R 2 ) 
and l'-2'-3' (for R 3 oR 2 ) . Explicitly, we obtain : 



If now we identify the ends of the interval ,i.e.O and 1 ,then both Eq.(3.1) 
and (3.2) are translations by a along the circle. The theory of circle home- 
omorphisms is well documented [17, 18] and, for the sake of continuity, we 
would like to remind to our readers several important facts which will help 
us to clarify the results just obtained and those which will follow. 

Let us notice that the circle S 1 is the quotient : S 1 = R/Z , i.e. the circle 
is the factor group of real numbers R modulo integers Z .Let / : S 1 ^S 1 
be some circle homeomorphism and let F : R — > R be a lift of / (e.g. 
/(S 1 ) = exp{27iia;} and F : x' = x + a ,that is F are translations in R).The 
rotation number t(F) can be defined via equation : 



[RioRaKx) = a - (1 - x) , 



(3.1) 



[R 3 oR 2 ](i) =a + l-(l-x) 



(3.2) 



t(F) = lim -(F n (x) 

InHoo n 



x) . 



(3.3) 



To define F(x) we notice that ,in general, 



F(x) = x + a(x), 



a(x + 1) = a(x), 



< x < 1 



(3.4) 



In terms of a(x) the rotation number can be written as well as 




Let the mapping F has a fixed point x* for some q: 



F q (x*) =x* + a q (x*) . 



(3.6) 
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It can be shown that, 



a q (x*) =p 



(3.7) 



and 



t(F) 



P 



(3.8) 



That is q translations in the lifting space are equivalent to "pure" rotation of 
the circle by the angle d = 27r| .Accordingly ,if Eq.(3.6) does not have fixed 
point (s) for any q,then r(F) = a where a is irrational. With respect to the 
situatioin depicted in Fig. 13 we have to conclude ,that the "motion" along 
the line l-2-3,etc will never stop or cross another line, i.e. we have obtained 
the foliation which is actually a labyrinth. Let us take now into account 
the P-H theorem and make the situation more realistic. To this purpose we 
attach the "tail" to S 1 semicircle and gently open the foliation leaf starting 
at a in Fig. 13. The result of such opening procedure is depicted in Fig. 14 
which represents the simplest example of a labyrinth: all incoming leaves are 
being trapped inside the labyrinth forever. 

The picture just described can be wasty generalized with help of the 
following auxiliary observation. Consider some nonintersecting line £ (closed 
or not) and place on the top of it a finite collection of foliated half cir- 
cles(semicircles) provided that none of these half circles are touching each 
other (unless the otherwise is specified). Let us place the remaining half discs 
on the bottom of £ in the way depicted in Fig. 15. We have obtained in this 
way an example of a "meander" labyrinth. Let us explain now what actually 
the word "meander" means. According to Ref. [34] ,a meander of order n is a 
closed nonselfintersecting curve which intersects another straight line in ex- 
actly 2n preassigned points(more exactlyan equivalence class of such closed 
curves which leave the straight line £ fixed). Evidently ,the line £ need not 
be straight in general.A few examples of the meanders of lower order and 
the corresponding meandritic labyrinths are depicted in Fig. 16. Although the 
four last meandritic labyrinths are degenerate (two half discs are touching 
each other), this may be still permissible (if other discs are not touching each 
other), e.g. see Example 2 of Ref. [33] .This is so because we still have to attach 
"tails" to some of the discs (e.g. see Fig. 14) to make them Y-type in order 
to be in accord with P-H theorem. Trough the tails the external lines pene- 
trate into the labyrinth and become trapped forever. Moreover, because of the 



Fig. 14 



Fig. 15 



Fig. 16 
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tails, the symmetry between different labyrinths depicted in Fig. 16 becomes 
broken so they all have to be consideredjust like the underlying meanders. 

The meandritic number M ra (that is the total number of meanders of or- 
der n) provides the degeneracy factor for the partition function of defects to 
be discussed in section 5. Its value is essential in the description of the di- 
clination binding-unbinding "melting" transition [14]in the non -Coulombic 
two dimensional "plasma" of such defects. Actuallythe above description 
is still incomplete since ,as it is well known from statistical mechanics of 
two dimensional plasma[21] ,the clustering effects should also be taken into 
account. This ,in particular,means that,instead of just one connected mean- 
der, we may as well consider a collection of non-intersecting meanders. Specif- 
icallyfor the multicomponent meander of order n we may choose k noninter- 
secting lines (l<k<n) so that the total number of crossings on all lines is 
2n. Accordinglyone can introduce the meandritic number to account 
for this fact. [35, 36] .This number will be estimated and used in section 5. In 
the meantime, it is useful to introduce ,yet another, set of meanders, the pro- 
jective meanders [19,37] .These are defined as follows. Consider 2n points 
on the circle S 1 . Divide the points into n pairs so that the chords connecting 
points in each pair will not intersect. Identify the diametrically opposite points 
on the circle,thus turning the disc into the projective plane(or sphere). Such 
formed set of curves in the disc is called the projective meander of order 
n.Evidentlythe number of possible projective meanders of order n is equal 
to the n-th Catalan number 



° n ~ (» + l)!n! (3 ' 9) 

which obviously reflects the combinatorics of the problem. We shall need the 
notion of the projective meanders in order to demonstrate that the above 
labyrinth construction is not artificial but , actually ,is intrinsic for the whole 
description of the measured foliations on surfaces. In addition,there is yet 
another way to look at the whole problem for the case of line fields. Because 
of its potential usefulness for the problems which involve 2+1 gravity ,we 
discuss this other approach in the next section. 
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4 Train tracks and pseudo- Anosov homeomor- 
phisms 

The logical development of dynamical systems on surfaces goes from the 
detailed treatment of the circle maps through the torus maps and then, to 
consideration of flows (foliations)on the Riemann surfaces of genus higher 
than one .Nevertheless, already the circle maps exhibit all qualitative fea- 
tures of more complicated situations. Moreover, as it was shown by Thurston 
[38] ,all surface homeomorphisms can be actually related to the circle maps 
through the construction which he calls the "earthquake". 

Let us briefly explain why this is so. (In the Appendix of Part II these 
arguments will be extended to 3 dimensional hyperbolic manifolds). As in the 
case of the circle maps discussed in the previous section,the surface homeo- 
morphism / : R — > R' (where both R and R' are some Riemann surfaces) 
can be lifted to the universal covering space which is either the Poincare 
upper halfplane H defined by 

H={z=x+iy G C | y > 0} (4.1) 

or the open unit disc D 2 

D 2 = {w=u+iv G C | u 2 + v 2 < 1} (4.2) 
which is related to H via mapping 

w=— , z G H . (4.3) 
z+i v 1 

Since every Riemann surface R of genus g greater than one can be obtained 
as the quotient R =H/T where V is some Fuchsian group, i.e. the group of 
linear fractional transformations of the type 

7 ( z ) = ,ad-bc = l (4.4) 

with a, b, c, d being real numbers ,we can always lift the flow on surface 
R to the flow on D 2 where by means of earhquakes the description for flows 
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becomes connected with the description associated with maps of the circle. Let 
now C be some closed curve on R which is homotopically nontrivial(that is 
it cannot be shrunk to the point). It can be proven 39 that, 

Theorem 4.1. C is freely homotopic to a unique closed geodesic I . 

Also, it is known [40] that, 

Theorem 4.2. For a closed Riemann surface of genus g there are 3g-3 
independent closed geodesies (please, see Fig. 5 of Part II ). 

Finally,following Thurston [25] and Casson and Bleiler [41] ,we introduce 

Definition 4.1. Disjoint union of geodesies is called lamination( L) ;the 
geodesies contained in L are called leaves of L . 

It can be very easily shown [17, 18] that, when lifted to H, the geodesies 
look like that depicted in Fig. 17. That is they are either halfcircles with center 
on the real axis or the lines perpendicular to the real axis .When transformed Fig. 17 
from H to D 2 ,the halfcircles go to halfcircles and the vertical lines to the 
diameters of D 2 as depicted in Fig. 18 .If we consider the case of lamination Fig. 18 
C then,by construction,the geodesies are nonintersecting on R and, whence 
on D 2 [42]. This means that in D 2 they look like the projective mean- 
ders(before identification of points on S 1 ) described in the previous section. A 
simple minded permutation of these geodesies in D 2 produces the Catalan 
number C n ,Eq.(3.9). In general, the situation is considerably more compli- 
cated [25,41] . 

If for a given Riemann surface of genus g with n punctures we have the 
inequality 2g+n-2 >0,then the mapping class group M(R) is defined through 
its action on £ .If / e M(R) and C = {ci, ...,c r } G C ,then the mapping / 
is reducible if 



Alternativelyit is being said that / is reduced by C. 

If, in accord with the theory of braids [43] ,we assume that Eq.(4.5) holds 
for any permutation of the set C (and this is indeed the case [41]) ,then the 
Catalan number ,Eq.(3.9) , naturally emerges The reducible case is not the 
only possibility however. If, say, / is reduced by C ,then let S = {si ,...,s m } 
be the set of components of the complement of R — C e.g. see Fig. 5 of Part 
II . Since / permutes Si ,let Hi be the smallest power of / such that 



f(C) 



c 



(4.5) 



(4.6) 
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For any rii > such mapping is called periodic .It could as well be that 
Eq.(4.6) does not have fixed point(s) for any finite rij. In this case, if Eq.(4.5) 
still holds, then the mapping is still reducible , but if neither (4.5) or (4.6) hold 
, then the mapping is called pseudo-Anosov. This term will be explained 
below. Before doing so, we notice that Eq.(4.5) can be lifted to D 2 so that with 
respect to the ends of geodesies lying on S 1 we have some sort of a mapping 
of a circle analogous to that discussed in the previous section. According to 
Thurston [38] ," every homeomorphism of the circle S 1 extends to the home- 
omorphism of the disc D 2 ".(In the Appendix to Part II we shall discuss the 
extension of this result to the case of 3-manifolds). Evidently, the reverse 
should be true as well. In which case , the lift of Eq.(4.6) to D 2 has its ana- 
logue in Eq.(3.6) and the case when Eq.(3.6) does not have solution for any 
q corresponds to the pseudo -Anosov type of mapping. Let us now explain 
the meaning of the word "pseudo-Anosov". According to Refs. [17,18] ,the 
following theorem can be proven. 

Theorem 4.3. The geodesic flow on H/Y is the Anosov flow. 

The best way to explain this is again pictorial .We lift the geodesies from 
R to H in order to obtain the picture given in Fig. 19. Accordingly, we can also Fig. 19 
look at the same picture by using the disc D 2 model. In this case we have the 
situation depicted in Fig. 20. The circles tangent to x-axis at some point x (in Fig.20 
H-plane) are called horocycles.lt can be shown, that the set of horocycles 
(touching x-axis at x ) and the set of geodesies (emanating from x ) are 
mutually orthogonal (thus forming stable and unstable manifolds, see below). 
For the geodesies directed as in Fig. 19 it is intuitively clear that the flow is 
expanding(unstable) .If we would change the direction of geodesics(i.e.if we 
change the time evolution from t to -t),then the flow will be contracting.lt 
can be shown 18 that the rate of expansion is e while the rate of contraction 
is e _1 (e=exp).Now ,we are ready for a general definition of the Anosov flow. 

Definition 4.2. If A : M-^M is some diffeomorphism of compact manifold 
M such that the tangent bundle T X M for any point x GM is decomposable 
into direct sum © of the type 

T.M^eY, (4.7) 

so that in some (Riemannian) metric ,and for some number A > 1 
the following inequalities hold 

\\A£\\ > A||e|| V^GX, and \\A^\\ < \' 1 \\ V \\ V^gY, (4.8) 
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where A* is the differential of the operator A (acting in tangent space) and 
\\...\\ is the usual square of the length in some (Riemannian) metric space, then 
the flow T X M is Anosov. 

Remark 4.1. If ,in addition, the flow contains some singularities, then it 
is called pseudo- Anosov[AA\. The sets X x and Y x are called (un) stable folia- 
tions respectively. 

Remark 4.2. It can be shown,that the geodesic flow on H is ergodic[18] 
. That is there is some measure ji{x) such that 



where A 1 is some evolution operator (continuous analog of A* in (4-8)). 

In the case of pseudo- Anosov flows it is customary to introduce the trans- 
verse (vertical) measure which can be intuitively understood using the fol- 
lowing picture, e.g.see Fig. 21 .That is, if locally the foliation looks like the Fig 
set of horizontal lines(away from singularities), Fig. 21a), then, when we try to 
restore the global picture,the local charts should be glued in such a way that 
the vertical distance A between the leaves on one chart is in agreement with 
the vertical distance A on another. If this is the case, the foliation is called 
measured .We shall be dealing only with measured foliations [10, 25] in the 
rest of this work in accord with Poenaru[10] and Langevin [11] . 

The above nice picture leaves us with no clues about the mechanisms 
by which the flows can be changed from the reducible to periodic or from 
the reducible to pseudo- Anosov,etc. As in the case of circle homeomorphisms 
(section 3), there must be some physical reasons of changing from one regime 
to another. The best way to get some physical feeling of the processes which 
cause "phase transitions" (i.e. changes in the flow regime) is through the no- 
tion of "train tracks" introduced by Thurston [20] and subsequently de- 
veloped by many other mathematicians,e.g.see Ref.[45].The idea behind the 
train track is simple but very powerful. To appreciate it , it is helpful to 
recognize that there is some similarity between the way knots(or links) are 
described and the train tracks. In the case of knots (links) one usually 
studies a "shadow" of a knot obtained by projecting it onto some arbitrary 
plane[46].The projection is a four-valent planar graph K without the dangling 
("dead") ends. Evidently, different knots may have the same projection for a 
plane which orientation is fixed. To distinguish between different knots(links) 



r 




(4.9) 
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one has to resolve each 4-valent vertex (i.e. to decide which strand is "over" 
and which is "under") and, in addition,to use the set of isotopy moves (the 
Reidemeister moves) to bring one projection in accord with the other. It is 
believed, that ,at least in principle,one can decide if two knots are equivalent 
by performing some finite (but could be very large !) sequence of the Reide- 
meister moves. In all these operations the physical nature of a knot is not 
playing any role. In the case of train tracks one is also dealing with a graph T 
without dead ends and also there are the isotopy moves(analogous to the Rei- 
demeister moves), but, in addition, and this is the most important, there 
are moves which do not respect isotopy. These moves change the topology 
of the graph T and could be associated with some physical processes as we 
shall explain below and in Part II (section 6 and Appendix). The graph T 
has two basic building blocks which can be easily recognized,e.g. see Fig. 22 Fig.22 
.One can see from this figure that, instead of having all leaves and precise 
angles ,one can actually "survive" only with topologically equivalent objects 
which are smooth at joints. With the basic building blocks just defined ,we 
can construct our first train track. In particular, instead of having rather com- 
plicated foliation pattern depicted in Figs. 13, 14, we can draw the following 
graph depicted in Fig.23. Of course, the r.h.s.of Fig. 23 provides still only a Fig.23 
coarse grained foliation pattern thus creating an illusion that the line is closed 
(and this is indeed the case for the "rational" labyrinth). The numbers on the 
graph are quite arbitrary but are subject to one restriction: at every switch 
the sum of "entering" numbers is equal to the sum of the "exiting " numbers 
.The numbers are associated with the invariant transverse measures which 
were formally introduced earlier. Since one can construct from the collection 
of 3-valent vertices the vertex of more complicated nature ,e.g. 4-vertex ,5- 
vertex,etc, it makes sense to talk about given vertex v& (k=l-M) in general.. 
If {af 1 } is the set of input branches and {a° ut } is the set of output branches 
,then for a given vertex v& the switch condition can be written as: 



The switch condition is analogous to the Kirkhoff rule for currents known 
in physics [47] .It could be proven [20, 44] ,that the assignment of weights 
(transverse measures) /i(a«) for each branch of the train track which is sub- 
ject to the switch conditions ,Eq.(4.10), at each vertex is equivalent to the 
reconstruction of the entire measured foliation (up to isotopies and White- 




(4.10) 
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head moves, Fig. 4) [23]. The new feature which makes the train tracks more 
complicated than knots (links) lies in additional moves which, unlike the 
isotopy,do change the topological type of the track.These moves have the 
major physical significance as it will be demonstrated shortly. Before doing 
so ,we would like to provide a list of topology changing moves .They are 
depicted in Fig. 24. The shift operation can be performed without reference Fig 
to weights but,of course, the switch condition, Eq. (4. 10), should be obeyed be- 
fore and after the shift. The split (or collapse) is dependent , however, upon the 
particular distribution of weights. The shift does not destroy (or create) the 
existing singularities while the split (collapse) may destroy (create) the sin- 
gularities thus apparently violating the P-H theorem. This deficiency can be 
easily corrected if some physics is taken into consideration .To do so, several 
steps are still required. First, following Ref. [45] ,let us consider the process 
of collapse as the matrix operation. That is, we form the column vector 
made of weights a', b', e' and look for the transition matrix S _1 such that 
x = S _1 x' or, more explicitlywhen d < a' ,we obtain, 



b 

c 
d 



/ 1 \ 

10 

10 

1 

\1 1 1/ 



/ a'\ 
b' 

c' 
d! 

W J 



(4.11) 



This result can be easily understood based on Eq. (4.10) and Fig. 24. 
Indeed, according to Eq. (4.10) we have a'+e' = c' (d > a') and also e'+d' = b' 
for the split diagram while a + b = e and c + d = e for the collapsed diagram. 



At the same time, using Eq.(4.11), we obtain: 



a ' ; b = b';c = c';d = d! 



but e = a' + d' + e'.The last result becomes an identity if the Kirkhoff sum 
rules Just obtained, are utilized. Evidentlythe corresponding matrices can be 
constructed for all elementary processes depicted in Fig.24 and,furthermore,if 
,for example, S is the split matrix, then S -1 is the corresponding collapse 
matrix,etc. Second, instead of considering the individual (local) processes 
,we can consider the state vector of the entire graph T which is just a column 
vector of all transverse measures which we denote as fl .The evolution of the 
entire graph is determined then by some transition matrix T so that 
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(4.12) 



Third ,this equation is still unphysical.To make it physical,we have to ascribe 
certain statistical weights to each of the diagrams depicted in Fig.24. The 
weights could be chosen either on the basis of statistical mechanics arguments 
(e.g. in the case of liquid crystals)or on the basis of quantum mechanical 
arguments ( e.g. in the case of gravity). For instance, in the case of liquid 
crystals one can think of energy of the defect formation and destruction 
(for more details ,please,see the next section and sections 4,6 of Part II)as 
well as about the average energy (per defect) as a function of v (defined 
by Eq.(2. 5)). Irrespective of the specific form of the weights ,it is clear ,that 
the time evolution of the graph T can be described in terms of the master 
equation(as it is usually done in statistical mechanics) 



The matrix Wij is not symmetric unfortunately.This can be seen already 
from the Eq.(4.11).In the case of 2+1 gravity the above equation could ,in 
principle, provide an alternative to the Wheeler-DeWitt equation[48] .Eq.(4.13) 
may have time -independent solution which (in terms of discrete maps, e.g. 
see section 3) ,by analogy with Eq.(4.5) ,we may call reducible .It may as 
well have a periodic (in time )solution which we may (or may not) asso- 
ciate with the previously defined periodic(in space) case. Finally, it may 
not have any stationary or periodic solutions. This then will be indicative of 
the pseudo-Anosov type of evolution. The actual analysis of Eq.(4.13) would 
require a large scale numerical work and,evidently,some drastic approxima- 
tions will be required. Therefore, for the time being,we would like to explore 
yet another possibility. 

Before doing so, several remarks should be made. 

Remark 4.3. Since in the case of train tracks there are no dead ends,this 
concept should be applied to open surf aces, e.g. to the open disc D 2 ,with some 
caution, e.g. see Ref.[49]. 

Remark 4.4. Since the surface of non-negative Euler characteristic con- 
tains no train tracks[20,45] some care should be taken to by-pass this diffi- 
culty. This is discussed to some extent in the Appendix (see also Refs.[33] and 



Remark 4.5. Without use of the probabilities leading to the master Eq.(4-13), 
different outcomes of iteration had been studied already [45, 49, 50] .Clearly,the 




(4.13) 



[49]) 
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process is reducible periodic or pseudo-Anosov depending upon what sequence 
of moves, Fig. 24 ,had been used. For example, if only direct processes are being 
considered, such as the right (R) and the left (L) splits and the shift (S), then 
the matrix sequence could be associated with some words in the alphabet made 
up of 3 letters [50]so that ,in general,Eq.(4-13) describes a random walk on 
the mapping class group M(R). 

Incidentally ,the 3 letters alphabet is being used for the meanders [34]. 
Hence, we would like now to study how one can use the meanders in resolv- 
ing the types of dynamical behaviour. This is accomplished to some extent in 
the next section while in the Appendix we provide some evidence that the 
meandritic labyrinths are indeed associated with the pseudo-Anosov homeo- 
morphisms. 

5 Phase transitions in the system of mean- 
dritic labyrinths 

Phase transitions in the system of meanders was recently considered in Ref . [35] 
in connection with the problem of folding of polymers. At the same time,the 
arguments used in Ref. [35] could be traced back to the arguments made by 
Peierls in connection with phase transintion in 2 dimensional Ising model 
[51] .To facilitate reader's understanding of the development which follows 
,we would like to provide here a summary of the arguments made by Peierls. 

In the case of planar Ising model one is expecting to have droplets of spins 
of the "wrong" sign embedded among the "bulk" which has the "correct" 
sign (e.g. "+" for the correct and "-" for the wrong). The interface between 
the correct and the wrong spins determines the interfacial energy Ex, of the 
droplet so that for the droplet of perimeter L the energy is given by El = LJ 
with J being known ,in principle, constant. The partition function Z/ for the 
Ising model can now be written as follows, 

00 

Zj = 2 G L exp{-pjL} , /T 1 = k B T, (5.1) 

L=0 

where ksT is the usual temperature factor. Since the Ising model is defined on 
the lattice ,L is necessarily discrete and Gl is the combinatorial entropic fac- 
tor which determines the number of ways the droplets of the total perimeter 
length L can be arranged on the lattice. 
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The difficulty in calculating Z/ lies in determining the factor Gl correctly 
.If, following Peierls,we make a drastic approximation :Gl ~ 4 L ( since for 
the square lattice the coordination number is 4), then Zj can be calculated at 
once with the result 



, where A = exp{— f3J}. The above expression makes sense only for A < 0.25 
so that 0.25 determines the critical temperature of the Ising model through 
equation : 4A* = 1 .This result is interesting to compare against the exact re- 
sult [52] :A* = \J >/2 — 1 ~ 0.64 (for J = 1) . Obviouslythe value 0.25 appears 
to be too low even if compared with the Bethe approximation [53] which yields 
3 (A*)' = 1 .This was noticed already by Peierls 51 .Nevertheless,the above 
crude estimate can be systematically improved and the Peierls arguments 
are the most powerful tool in general study of phase transitions in discrete 
spin systems as is well known [53] .The exact result ,(A*) 2 = \/2 — 1, is valid 
only for the square lattice. In the case of triangular lattice the exact result is 
different: (A*) 2 = (x/S) , while for the hexagonal lattice ,it is different again : 
(A*) 2 = (2 + a/3)" 1 , e.g. see Ref.[52].Hence,the critical temperature is actually 
a function of two variables :J and z ,where z is the coordination number of 
the lattice. Something similar will happen in our (meandritic )case too. 

In the case of Ising model we can effectively define the notion of an or- 
der parameter. This is less obvious in the case which we are going to con- 
sider. What is important for us, however, is the fact that the partition func- 
tion Z/ diverges. We shall ,by analogy with the Ising model , associate such 
divergence with the phase transition, in our case from the meandritic(pseudo- 
Anosov, according to Ref. [33] and Appendix) to the periodic and/or reduced 
phases (since these are the only existing possibilities). 

Using the results of section 2 ,Eq.(2.5) and section 4 of Part II, Eq. (4. 49), the 
total surface energy of n defects can be estimated as 



where the charge q is equal to the index of the singularity, Eq. (1.11), and K 
is some known constant (defined in section 4 , Part II) related to the surface 



Zi ~ (1 - 4A) 



-i 



(5.2) 




(5.3) 
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tension. In arriving at Eq.(5.3) we had assumed that only Y's and thorns are 
present (so that the magnitude of charges are all the same). To estimate the 
combinatorial factor (analogous to Gl in Eq.(5.1)) the following arguments 
are helpful. 

According to Ref.[36] ,any meander can be built by superposition of two 
arc configurations of the same order :one is considered to be the top while an- 
other is the bottom as depicted in Fig. 25. Since both the top and the bottom Fig.25. 
are configurations of the same order, then by concentrating our attention,say 
,on the top configuration, we can obtain a projective meander (e.g. see sec- 
tion 3 )by identifying the beginning and the end of the horizontal straight 
line. In this case we already know that the number of possible configurations 
is C n .Evidentlywe can do the same for the bottom arc system thus obtain- 
ing another C n . If we try all possible permutations, then, a priorythere is no 
guarantee that if the top and the bottom are connected together we shall 
obtain just one connected meander. More likely ,we shall obtain C\ meanders 
some of them connected and some not. More exactly, let M(x) be the meander 
generating function, i.e. 

oo 

M(x)=^M n x n , (5.4) 

n=0 

where M n was defined in section 3. There is the following 
Theorem 5.1. C n <M n < C 2 n 

Proof .The upper bound was just established. The lower bound can be 
established if we associate the system of meanders with the 3 letter alpha- 
bet (as was briefly mentioned in section 4). For more details,please, consult 
Ref.[34]. □ 

Among all of C\ meanders of order n not all are topologically different. Let 
N n denote the number of topologically distinct configurations of meanders 
which altogether pass through 2n points (which actually may belong to dif- 
ferent lines), then the generating function for these numbers could be defined 

as 

oo 

N(x)=^N n x n . (5.5) 

n=0 

If we disrespect the topological differences, then the generating function for 
the system of meanders can be written as 
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oo 

B(x)=£c2x». 

?1=0 



(5.6) 



Theorem 5.2. Functions B(x) and N(x) are connected with each other 
through the functional equation 

B(x)=N(xB 2 (x)) (5.7) 
Proof.Please, consult Ref.[34]. □ 

(k) 

Corollary 5.2 . The meandritic numbers M n introduced in section 3 are 
generated via the following generating function [35, 36] 

n 

m„(g)=^Mi fc )g fc (5.8) 

fc=0 

provided that m n (g=l)=C^ l . 

Based on the results just described, we are ready now to write down the 
partition function for the system of meandritic labyrinths. It is given by 

oo n 

W = E E M n^ kxn > x=exp{-/3J}, (5.9) 

n=0 k=0 

where J is defined by Eq. (5. 3). This equation ,in principle, is exact but, it re- 
quires knowlege of g. Evidently^ is the fugacity which determines the average 
number of meanders in the cluster <k> which is given by 

<k>=gj|ln^(x) . (5.10) 

This quantity is hard to estimate (see,however ,some attempts in this di- 
rection below). And, even if we would succeed, still,we would need to invert 
the infinite power series in order to write the fugacity g in terms of <k> . 
Therefore,we would like to make an approximation based on the exact result 
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coming from Eq. (5.8) for g=l.Thus ,we obtain the following approximation 
for the partition function of the system of meanders: 



oo 

Z(x)=J2cy , x=exp{-/?J}. (5.11) 

n=0 

Remarkably enough, this partition function admits the exact resummation[34] 
with the result (t 2 =x) : 



2tt 



Z ^ = 4^2 + 27 / dVy/ 1 ~ 8t cos V 9 + 16t ' j • ( 5 - 12 ) 

At the same time,one can study the convergence of Z(x) based on known 
asymptotic value for C n : 
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C n ~ const— , n — > oo. (5.13) 



From here,we obtain, 



% ~ 16 (5.14) 

and,therefore,0<x< ^ .Substitution of x* = 16 into Eq.(5.12) produces fi- 
nite result : 



For x> ^ the partition function diverges, of course. 

Clearlythe condition 16x* = 1 determines the critical temperature for the 
fixed value of v or determines the critical "density" v* for a fixed temperature 
and the value of surface tension K in view of Eq.(5.3). 

The analysis presented above is still incomplete since one can calculate 
now N(x) in view of Eq.(5.7). Following Ref.[34] ,let us multiply both sides 
of eq.(5.7) by t (t 2 =x).We obtain, 
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$(t)=tB(t 2 ) = tN($ 2 (t)) . (5.16) 
Let, furthermore, A = <3>(t) so that t=$ _1 (A). Then,we obtain, 



N(A2) = ¥=W ' (5 ' 17) 

From here it follows, that since using closed form result , Eq.(5.12),<l>(t)=tZ(t 2 ) 
can be calculated, at least in principle, then N(A 2 ) can be calculated as 
well.The radius of convergence for N(x) is equal to (^ 2L ) 2 .Obtained results 
are in formal qualitative accord with that known in physical literature on 
liquid crystals [14]. But now we know that the disordered (liquid-like) phase 
is actually pseudo-Anosov while the ordered (hexatic) is periodic and the re- 
ducible is solid-like. The phase transition mechanism , however, has absolutely 
nothing to do with the Kosterlitz-Thouless type of transition which was dis- 
cussed in section 2. Additional details related to this subject are provided in 
section 4 ,Part II. 

Note added in proof: When this paper was completed we had found a 
paper by Ishikawa and Lawrentovich (Europhys. Lett. 41,pp. 171-176 (1998)) 
in which the Whitehead moves depicted in our Fig.4 had been detected ex- 
perimentally. In addition, the paper by Penner (Adv. in Math.l01,pp31-49 
(1993)) had also came to our attention. In this paper the train-tracks are be- 
ing used to describe folding of RNA. Since in Ref.[35] the folding of proteins 
was considered with help of the meanders, the above references provide addi- 
tional support to the results presented in sections 3-5 which relate meanders 
to train tracks. More details could also be found in section 7.5 of Part II. 

Appendix 

In this Appendix we would like to provide some arguments in favor of the 
statement made in section 4 that the meandritic labyrinths are associated 
with the pseudo-Anosov homeomorphisms. 

To this purpose let us consider the simplest labyrinth which is made out 
of two copies of the disc D 2 both of which containing just two thorns (to be 
in accord with P-H theorem). When glued properlythese two discs will form 
a sphere S 2 with foliation forming a labyrinth (e.g. see Ref.[33],page 29). We 
need now to explain why such foliation could be of pseudo-Anosov type. 
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Consider an auxiliary problem about the foliations on the torus[17, 18, 41] 
.If we regard the torus T 2 as a quotient of R 2 by the integer lattice Z 2 ,then 
the homeomorphisms h a of T 2 are generated by the group GL 2 (Z) since any 
element a of GL 2 (Z) maps Z 2 into itself thus inducing the continuous map 
h a : T 2 ^T 2 .The homeomorphism is orientation preserving if det (a) = 1 
(Vet eGL 2 (Z)) and, in this case ,GL 2 (Z) becomes SL 2 (Z) which is represented 
by 2x2 matrix given by 



The eigenvalues of the above matrix are obtained, as usual, as the roots of 
the characteristic polynomial 



There are three possibilities in general: 

1) the roots of (A. 2) are complex (when a+d=0,l or -1); 

2) both roots of (A. 2) are equal to ±1 (when a+d=2 ); 

3) the roots of (A. 2) are distinct reals (when \a + d\ > 2) . 

By analogy with the circle maps considered in section 3, it can be 
shown [16, 41], that the first possibility produces periodic maps of T 2 ,the 
second produces reducible maps of T 2 (that leave a simple Jordan curve 
fixed (possibly with reversed orientation)) while the third possibility is re- 
sponsible for the Anosov (not pseudo-Anosov) type of flow [17, 18] on T 2 
which was defined in section 4. Consider now our S 2 ,made up of two copies of 
D 2 ,each containing two thorns. According to the results of Part II (section 
3) thorns are singularities in the complex z-plane (or on S 2 ) of the type 
{\fz) 1 .As it is usually done in the complex analysis, we can make branch 
cuts (as depicted in Fig.Al)in order to construct the two sheeted Riemann Fig.A.l 
surface. In Fig.A.l these two sheets are depicted as onion-like sphere with 
two layers. If the outer sphere is peeled off and placed opposite to the inner 
sphere, as depicted in Fig. A. 2, then,by obvious identification of edges ,the sur- 
face which is homeomorphic to T 2 is formed. Hence, T 2 is a two-fold ramified 




(A.l) 



t 2 - (a + d)t + (ad - cb) = 



or, in view of (A.l) , 



t 2 - (a + d)t + 1 = 0. 



(A.2) 
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branched covering of S 2 .The points of ramification are the above thorn-like 
singularities. In general,one can prove the following 

Theorem A.l. Every closed orientable surface (of finite genus) is two- 
fold ramified branched cover of S 2 . 

Proof .Please, consult Refs. [55-57]. □ . 

Finally, if for T 2 we have the case of the Anosov flow, then it is being 
transferred onto S 2 in the form of pseudo-Anosov flow (because of 4 singu- 
larities placed on S 2 ).More exactly,the foliation on S 2 is the quotient of the 
Anosov homeo(or diffeo)morphism on T 2 via 2-fold ramified branched cov- 
ering of S 2 by T 2 . For more than 4 defects in view of Theorems A.l and 
4.3 of the main text ,it is also possible ,in principle, to establish connection 
between the labyrinths and pseudo-Anosov homeomorphisms. Please, consult 
Expose/11 and Expose/12 of Ref[23] and section 4 of Ref.[44] for additional 
details. For more information on branched coverings,please, consult Ref. [57]. 
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Figure Captions 



Fig.l Typical singularities of the vector fields in the plane 

Fig. 2 Typical singularities of the line field 

Fig. 3 The field lines around the dipole a) and the source b) 

(In the case of a sink the directions of lines are changed into the opposite) 

Fig. 4 The allowed Whitehead moves for the line fields(" creation") 

Fig. 5 The additional configurations leading to " destruction" of defects 

Fig. 6 Possible "periodic" phase in the case of minimal number of line field 

defects in the disc 

Fig. 7 Possible "periodic "phase made up of 14 defects in the disc. 
Fig. 8 Possible hexatic-type phase obtained in the thermodynamic limit 
Fig. 9 The elementary triangle used in P-H theorem calculations related 
to the hexatic-type phase depicted in Fig.8 

Fig. 10 a) Another possible hexatic phase which is permitted by the Euler 
theorem; 

b)The elementary triangle used in Hopf-like calculations of the 
Euler characteristic 

Fig. 11 Two dimensional plane can be covered by hexagons. At the same 
time, one can cover the plane by the triangular lattice so that for the 
complete coverage 3 hexagons should fit together around each vertex 
Fig. 12 Simplest "building blocks" for measured foliations 
Fig. 13 Simplest nontrivial global foliation pattern made of simple 
building blocks 

Fig. 14 The structure of disordered phase (labyrinth) in the case 
of minimal number of line defects (to be compared with the 
ordered phase depicted in Fig. 6) 

Fig. 15 The simplest example of the meander labyrinth 
Fig. 16 The simplest meanders of order n=l-3 
Fig. 17 Geodesies on the Poincare upper H- plane 
Fig. 18 Geodesies on the unit open disc D 2 
Fig. 19 Anosov flow on H-plane 
Fig. 20 Anosov flow on the disc D 2 

Fig. 21 For the measured foliation the vertical distance between the 
leaves on two adjacent patches (local charts) must be the same 
Fig. 22 Basic building blocks of the train track 
Fig. 23 A simple train track and the foliation pattern which is 
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associated with it 

Fig. 24 Major topology nonpreserving moves for the train tracks 

Fig. 25 Construction of a typical meander 

Fig.A.l Two sheeted covering of S 2 before glueing 

Fig. A. 2 Identification of cuts between the inner and the outer spheres 

produces torus T 2 
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